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In HISQ simulations by the MILC and Fermilab Lattice collaborations, both the light quarks and 
the charm quark are staggered. We extend staggered chiral perturbation theory (SZPT) to include 
Q_i such all-staggered heavy-light mesons. We assume that the heavy quark action is sufficiently 

improved that we may take aniQ « 1 (where ihq is the heavy quark mass), but also that niQ » 
Agco so that a continuum heavy quark expansion is appropriate. Using this SZPT, the leptonic 
y—i decay constant of the heavy-light meson is calculated at next-to-leading-order. The pattern of 

taste splittings in the heavy-light meson masses is also investigated. 
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1. Introduction 

Heavy-light meson systems provide some of the best ways to test the Standard Model and 
look for signs of new physics. Lattice QCD provides a means of carrying out non-perturbative 
calculations from first principles and with controlled errors. In setting up a lattice QCD calculation, 
a key choice is the form of the lattice action for the quarks. 

Staggered fermions [1] are an efficient approach to simulating light quarks. The "highly im- 
proved staggered quark" (HISQ) action [2] makes it possible to treat charm quarks with the same 
action as the light quarks. Thus "all- staggered" simulations of D and D s mesons are now possible 
[3]. There are several advantages to this all-staggered approach. Probably the most important is 
that, since heavy and light quarks have the same action, there are partially conserved heavy-light 
axial and vector currents, which therefore need no renormalization. 

Lattice computations often involve an extrapolation in light quark masses to the physical up 
and down masses, and always require a continuum extrapolation in lattice spacing. A version of 
chiral perturbation theory (ZPT) that includes the effects of the discretization errors can help to 
control these extrapolations. Here, we develop chiral perturbation theory for all-staggered heavy- 
light mesons. 

Reference [4] works out a closely related chiral theory for heavy-light mesons with staggered 
light quarks but non-staggered heavy quarks (for example, Fermilab [5] or NRQCD [6] quarks). In 
that case, the doubler states of a heavy quark are treated as integrated out, and therefore, heavy- 
light mesons have a single taste degree of freedom associated with the light quark. Here, we need 
to extend the program developed in Ref. [4] to include staggered heavy quarks with a taste degree 
of freedom. We assume that the staggered action used (e.g., HISQ) is improved sufficiently that we 
can treat the heavy quark as "continuum-like," with small corrections from cutoff effects. We refer 
to this assumption in short-hand as taking aniQ <C 1, where itiq is mass of the heavy quark. We can 
then use the Symanzik effective theory (SET) [7] to describe the discretization effects on the heavy 
quarks, as well as on the light quarks. 

In the continuum limit, there is an exact SU (4) symmetry acting on tastes; this symmetry is 
broken at &{a 2 ) in the lattice spacing a. The corresponding discretization errors in the light-light 
sector split the masses of mesons with different tastes, which may be understood using staggered 
chiral perturbation theory (SZPT) [8, 9]. For typical values of a 2 , the taste splittings of the light 
pseudoscalar mesons can be comparable to the masses themselves. Schematically, we say a 2 ~ m|, 
where appropriate powers of Aqcd or A x (the chiral scale) are implicitly inserted to match the 
dimensions in such comparisons. The taste splittings are therefore included in the leading order 
(LO) light-light Lagrangian, which is of &(m~). 

For heavy-light mesons composed of staggered quarks, the situation is different. The LO La- 
grangian in the continuum is of 0(k), where k is the residual momentum of the heavy-light meson, 
and we assume that k ~ m n . Thus, it is reasonable to treat taste violations, which are of &{a r), 
as next-to-leading order (NLO) corrections, and that is what we do here. The LO heavy-light La- 
grangian is then taste invariant. This power counting can be checked with HISQ simulations, where 
the splittings in squared meson masses remain roughly constant as one increases a valence quark 
mass from the light quark regime to the charm regime [10]. (See also Fig. 1 below.) Therefore 
the splittings for the masses themselves are much smaller for heavy-light mesons than for light 
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mesons. For example, at a « 0.12 fm the measured [10] taste splitting between the root-mean- 
squared (RMS) D s meson and the lightest D s meson is only about 11 MeV, while it is about 110 
MeV for the pion (80% of m„). 

2. The SZPT Lagrangian for Heavy-Light Mesons 

We first discuss a heavy-light meson where the light quark has a taste degree of freedom, 
but the heavy quark does not [4]. Due to the heavy quark spin symmetry in the static limit, the 
heavy vector and pseudoscalar mesons are incorporated into the following field, which destroys a 
heavy-light meson, 

H a = l -^-[fBl a + iy 5 B a ] , (2.1) 

where v is the meson's velocity, and a is the combined light quark flavor-taste index. Note that 
we use B to denote a generic pseudoscalar heavy meson and B* to denote the corresponding vector 
meson, but the practical application of this calculation will, at least in the first instance, be to D 
mesons. 

Now, we assume the heavy quark is also implemented by a staggered fermion. Then we 
generalize the definition of the annihilation operator of a heavy-light meson as 

H aa = ^ [tBl aa + ijsBaa] , (2.2) 

where v is the meson's velocity, a is the heavy-quark taste index, and a is the combined flavor-taste 
index of the light quark. The conjugate field creates a heavy-light meson 

H aa = Y H} iaYo = [fB% a + iy&la] ^ . (2.3) 

So far H is treated as a 4 x An matrix in the taste and the flavor space of quarks. Instead of 
attaching separate indices to the tastes of the light and the heavy quarks, one might use just one 
index as the taste of the meson. Then, one treats H as a ^-component vector in the flavor space 
of the light quark, while each element (Hi) is a 4 x 4 matrix in the taste space of the meson. To 
implement this approach, the light flavor-taste index is first traded for a pair of indices representing 
flavor and taste separately. We use Latin indices in the middle of the alphabet as pure 

flavor indices. We can then combine the tastes of the heavy and the light quarks and use an index 
such as S (where E = 1, . . . 16) as the taste of the meson. Therefore, the /th element of the field 
destroying a heavy-light meson in the light flavor space can be represented by //, = £s=i Hj^Tz and 
its conjugate by //, = £gLi H^Tz, with the Hermitian taste generators 7V £ {^5, i%n V ,%n,%i}. 

The leading-order staggered version of the heavy-light meson Lagrangian then looks the same 
as in the case with only light staggered fields, or indeed in standard continuum XPT, with the only 
change being the extra taste degrees of freedom in the fields. The LO chiral Lagrangian involving 
the heavy-light meson fields is: 

JSfi = -iTr(HHv-D) + g7t Tr(HHY >1 r5^) , (2.4) 

where (HH) a b = H aa H a b, and Tr, here and below, means the complete trace over flavor and taste 
indices and, where relevant, Dirac indices. This Lagrangian has separate SU (4) taste symmetries 
on the heavy and light quarks, as well as spin symmetry of the heavy quark. 
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To derive the heavy-light decay constants, we also will need the chiral representative of the 
axial heavy-light current. Alternatively, one can work with the left-handed current. The left-handed 
current that destroys a heavy-light meson of taste S and light flavor i is which at LO takes the 
form 

Jto = | M 1 ^ (1 - K)HoW>) , (2.5) 

where tr^ is a trace over Dirac and taste indices, and is a constant row vector that fixes the 
flavor of the light quark: (X^)j = 8ij. The decay constant fs m is defined by the matrix element 

«is(v)) = ifBay/m&PSssiSa , (2-6) 

where the state |B,s(v)) is normalized non-relativistically, corresponding to our non-relativistically 
normalized heavy meson field B. At LO in the heavy-light chiral theory, Jlo" = itoPBtw, which 
gives /lo = K/y/m^. 

To go beyond leading order, we must encode the discretization errors in the chiral theory. 
To do this we consider a series of effective field theories. Assuming omiq « 1, we can derive 
a Symanzik effective theory (SET) from the staggered lattice Lagrangian. One can then use the 
fact that niQ is large compared to Aqcd to organize heavy quark effects with heavy quark effective 
theory (HQET). Finally, when residual momenta and light quark masses are small compared to the 
chiral scale A x ~ 1 GeV, the physics of light-light and heavy-light mesons may be described by 
a chiral effective theory. Here, we assume the power counting p 2 ~ ra 2 ~ m q ~ a 2 for the light 
mesons as in Ref. [4]. 

There are three different types of 4-quark operators in the SET at order a 2 : 

a 2 @ ss 'tt' = c\ a 2 qi(Y s ®%t)qiqi>(Ys<<S>%ti)qi> 
+ c 2 a 2 qi{%®%t)qiq h (Ys'®E>t>)qh 

+ c 3 a 2 q h (y s <g> £,t)qhq~h'{Ys' ® €f , ( 2 -7) 

where qi and q^ are the light and heavy quark fields, and y s and £ f are any of the 16 spin or 
taste matrices, respectively. The operators listed in Eq. (2.7) are generic, and stand for a linear 
combination of all operators with the same heavy and light quark structure. The first type of 
operator has only light quarks and hence respects the heavy taste symmetry. Its contributions to 
the NLO chiral Lagrangian are thus essentially the same as those that appear when the heavy quark 
has no taste, i.e., those of Ref. [4]. Although the third type breaks the heavy taste symmetry, it 
contributes only in trivial ways to the heavy-light Lagrangian, which by definition has only two 
meson fields and does not describe heavy-light meson scattering. New chiral terms arise only from 
from the second type of operator, which has both heavy and light fields. 

Using a spurion analysis, we find that the operators of the second type lead to the following 
two contributions to the mass term at NLO Lagrangian: 

J2g 2 = K A1 a 2 Tr(H^H^ 5 ) + K A2 a 2 Tr (H^H^) +K A3 a 2 Tr (H^H^ 5 ) 
+ K A4 a 2 Tr (H^ V H^) + K A5 a 2 Tr {Hy 5il Hf 5 ) 
+ K A6 a 2 Tr {Hy 5 ^ 5 Hy^^) +K A1 a 2 Tr (Ey^Hf* 1 ^) 
+ ^ A8 a 2 Tr(F7 MV ^ A //7 v ^ A5 )+ J fiC4 9 a 2 Tr(F7 5 ^ vA // 7 ' ,5 ^ v ) , (2.8) 
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and 



= I ^i« 2 Tr(// 7v/l ^/// v ^) +J fi: B2 a 2 Tr(//7 v ^5^ V ^5) 



+ K B3 a 2 vf i v^Tr(H^H^ v )+K B4 a 2 Tr(Hy 5 ^ Vfl HY^^ v ) , (2.9) 



where J^ 2 2 (the "type-A" contribution) results from operators that are invariant over the full Eu- 
clidean space-time rotation group, 50(4), as well as a corresponding 50(4) of taste, and 
(the "type-B" contribution) results from operators that couple spin and taste and break these 50(4) 
symmetries. 



3. Results and Conclusions 



The 6 '(a 2 ) contributions to the Lagrangian, J2^ 2 2 and Jz^a 2 ' §* ve different NLO mass correc- 
tions to different tastes of the heavy-light mesons. The type-A terms split the masses into the five 



50(4) taste multiplets, labeled by I, V, T, A, and P. The type-B terms split these representations and 
give different masses to the time and spatial components, such as £q an d %i for the vector taste (V) 
representation. Note that no splitting, either between 50(4) multiplets, or within the multiplets, 
comes from the one-loop diagrams. This can be proved using the exact SU (4) heavy-quark taste 
symmetry of the LO Lagrangian and the discrete taste symmetry, which corresponds to the shift 
symmetry in the staggered action [4] . 

Based on the pattern of splittings seen at LO in light mesons, the 4-quark operators with 
the taste structure ^5 and spin /, 75, or y^ v appear to be dominant. These operators contribute 
to the C4 chiral operator [9], which gives the characteristic P, A, T, V, I (lowest-to-highest) or- 
dering of squared masses, with roughly equal spacing. Note that, for light mesons, only type-A 
4-quark operators are relevant to LO, because type-B operators have no chiral representatives un- 
til NLO [8]. Here the corresponding type-A 4-quark operators give rise to the chiral operators 
KA3d 2 Tr(H ^SfiH^uy) and KA^a 2 Tr(Hy^ v ^ 5 xHj V11 ^5). These produce the same equal-spacing 
pattern for the 50(4) taste representations of heavy-light mesons: See Table 1. For type-B opera- 
tors, one might guess that the taste ^5 and spin y^ v 4-quark operator would be dominant, since it is 
the only type-B operator that has the same spin and taste as one of the dominant type-A operators. 
This 4-quark operator gives rise to the chiral operator ££2fl ,2 Tr(//y V|U <i;5^//y' JV § iU 5). Table 2 
shows the pattern of mass splitting that stems from this type-B operator. 

The patterns given in Tables 1 and 2 are qualitatively present in the MILC data, shown in 
Fig. 1. Note in particular the "sc" case, for which the errors are small enough that the pattern of 
50(4) breaking for heavy-light mesons is clear. It is non-trivial that the time component of taste 
is higher than the space component in two cases (§0 vs. and 5o vs. ^j) but not in the third case 
(§05 vs. ^5), just as in Table 2. Although the chiral theory is not applicable to the "cc" case, it 
is interesting to see that the structure that would correspond to the dominant type-B operator gets 
particularly strong there, with near degeneracies of between members of different 50(4) multiplets, 
in particular §0 and /, or §/o and 

We can also calculate the decay constant of the heavy-light meson to NLO. At one loop, we 
express the decay constant as 




(3.1) 
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Table 1: Taste splittings due to the apparently dominant type- A operators. 
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Table 2: Taste splittings due to the apparently dominant type-B operator. 
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Figure 1: MILC HISQ ensemble at a ~ 0.15 fm and m/ = 0.2w s [10]. Squared mass splitting between pions 
of different tastes and the Goldstone pion in units of r\ . The types of quarks in the mesons are shown on the 
abscissa: 1, s, and c stand for light (u,d), strange, and charm quarks, respectively. 



The lowest order term f\ depends only on the light valence flavor as /g = k/ \jrriB x because 
the taste splittings of heavy-light masses appear at NLO. We can divide the terms contributing to 
the decay constant at NLO into two parts. The first part comes from the operators which obey the 
SU (4) taste symmetry of the heavy quarks, while the second stems from those breaking the heavy 
taste symmetry. The first part contributes both to the analytic terms and to the chiral logarithms, 
5/g r in Eq. (3.1). It is independent of the taste of the meson, and indeed is the same as the result in 
Ref. [4], because it obeys both the residual discrete taste symmetry and the SU (4) taste symmetry 
of the heavy quarks at LO. The second part contributes only to the analytic terms in Eq. (3.1), 
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which depend on the taste of the meson. 

We find the decay constant of the heavy-light meson in the partially quenched case is: 



1 + 



1 1 + 3*2 



1 



167T 2 / 2 



16- 



m xf,E>) 



a 2 8y- £ 



Rf\^;^)l{m)) +[V^A] 



+ c s (m u + m d +m s ) + c v m x + c ai za , 



(3.2) 



where x is the valence flavor, S is the valence taste, / runs over the three sea quarks u, d, and s, S' 
runs over the 16 meson tastes, and other notation is explained in Ref. [4]. In Eq. (3.2), c a ,s is the 
only coefficient which depends on the taste of the heavy meson. 

In summary, we have generalized the chiral Lagrangian of a heavy-light meson to the case 
where both heavy and light quarks have taste degrees of freedom. We have obtained the NLO 
chiral Lagrangian, which breaks the SU (4) heavy taste symmetry. Moreover we have derived the 
NLO decay constants and the taste splittings of the heavy-light masses, and have used them to 
understand, qualitatively, the pattern of splittings seen in the heavy-light HISQ data. 

We thank our colleagues in MILC for Figure 1 , and Andreas Kronfeld for helpful discussions 
on the use of effective theories for this problem. Our work was supported in part by the U.S. 
Department of Energy under Grant DE-FG02-91ER40628. 



References 

[1] J. B. Kogut and L. Susskind, Phys. Rev. D 11, 395 (1975). T. Banks et al, Phys. Rev. D 15, 1 1 1 1 
(1977); L. Susskind, Phys. Rev. D 16, 3031 (1977). 

[2] E. Follana et al. [HPQCD Collaboration], Phys. Rev. D 75, 054502 (2007) [hep-lat/06 10092]. 

[3] E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu [HPQCD Collaboration and UKQCD 
Collaboration], Phys. Rev. Lett. 100, 062002 (2008) [arXiv: 0706. 1726]; J. Bailey et al. [MILC and 
Fermilab Lattice Collaborations] , P o S ( LAT 2 012)159. 

[4] C. Aubin and C. Bernard, Phys. Rev. D73 014515 (2006) [arXiv:hep-lat/0510088v3]. 

[5] A. X. El-Khadra, A. S. Kronfeld and P. B. Mackenzie, Phys. Rev. D 55, 3933 (1997) 
[arXiv:hep-lat/9604004]. 

[6] G. P. Lepage and B. A. Thacker, Nucl. Phys. Proc. Suppl. 4, 199 (1988); B. A. Thacker and 
G. P. Lepage, Phys. Rev. D 43, 196 (1991). 

[7] K. Symanzik, Nucl. Phys. B 226, 187 (187). 

[8] W. Lee and S. Shaipe Phys. Rev. D60, 114503 (1999) [arXiv:hep-lat/9905023]. 

[9] C. Aubin and C. Bernard, Phys. Rev. D 68, 034014 (2003) [arXiv:hep-lat/0304014] and Phys. Rev. D 
68, 074011 (2003) [arXiv:hep-lat/0306026]. 

[10] A. Bazavov et al. [MILC Collaboration], in preparation. 



7 



